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ABSTRACT 



The stability of a dilute plasma to local convective and rotational 
disturbances is examined. A subthermal magnetic field and finite thermal 
conductivity along the field lines are included in the analysis. Stability 
criteria similar in form to the classical H0iland inequalities are found, but 
with angular velocity gradients replacing angular momentum gradients, and 
temperature gradients replacing entropy gradients. These criteria are indifferent 
to the properties of the magnetic field and to the magnitude of the thermal 
conductivity. Angular velocity gradients and temperature gradients are both 
free energy sources; it is not surprising that they are directly relevant to the 
stability of the gas. Magnetic fields and thermal conductivity provide the 
means by which these sources can be tapped. Previous studies have generally 
been based upon the classical ff0iland criteria, which are inappropriate for 
magnetized, dilute astrophysical plasmas. In sharp contrast to recent claims in 
the literature, the new stability criteria demonstrate that marginal flow stability 
is not a fundamental property of accreting plasmas thought to be associated 
with low luminosity X-ray sources. 

Subject headings: accretion, accretion disks — black hole physics — 
convection — hydrodynamics — instabilities — turbulence 
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1. Introduction 

Accretion onto compact objects is generally possible only if specific angular momentum 
is somehow extracted from fluid elements. It is now known that magnetic fields in a 
differentially rotating fluid cause a breakdown of laminar flow into turbulence (Balbus & 
Hawley 1998), and that this turbulence leads to a much enhanced angular momentum 
transport. Despite the fact that turbulence is to some extent amenable to detailed numerical 
simulation, the theory of turbulent transport lacks a good phenomenological description. It 
may well be that one does not exist. 

To make progress, some form of idealization is usually necessary, and a common 
approach has been the following: the magnetohydrodynamic (MHD) turbulent transport is 
modeled as an enhanced Navier-Stokes viscosity in an unmagnetized fluid. The idea is that 
the flow of this ansatz hydrodynamical fluid may then be subject to further instabilities. In 
particular, the enhanced viscosity is also an energy source, and heating in the inner regions 
of some types of accretion flow leads, it is argued, to a convectively-unstable temperature 
gradient. It has been further argued that, much like stellar convective zones, the flow would 
quickly evolve to a state of marginal stability. Proponents of these convection-dominated 
accretion flows ('CDAFs') suggest that this process explains the very low X-ray luminosities 
associated with many black- hole candidates (Narayan, Igumenschev, & Abramowicz 2000), 
because the inward angular momentum transport associated with thermal convection would 
precisely offset the usual outward MHD transport. In the stellar case, marginal stability 
amounts to adopting an adiabatic temperature proflle; in an accretion flow matters are 
more complex, with constant entropy and angular momentum surfaces satisfying marginal 
stability by the H0iland criterion (Quataert & Gruzinov 2000). 

We are therefore motivated to consider the equilibrium of a rotating, magnetized, hot 
dilute plasma. The results of this investigation may be surprising to the reader, because 
they are at odds with what has become a standard approach. Our findings are relevant to 
understanding current problems in X-ray accretion sources, but they are also of general 
fiuid dynamical interest. They are easily stated: (1) In a rotating, stratified, dilute plasma, 
the presence of any magnetic field and any Coulomb-based thermal conductivity renders 
the classical H0iland criteria insufficient for flow stability (to axisymmetric disturbances). 
Maximum growth rates are rapid (dynamical time scale), and insensitive to the fleld strength 
and the value of the conductivity. (2) Where entropy {S) and angular momentum (Z) 
gradients appear in the classical H0iland formulae, they must be replaced with temperature 
(T) and angular velocity {Vt) gradients. Spherical adiabatic accretion is unstable. (3) 
Convective instability is inseparable from rotational instability, in the sense that the 
breakdown of laminar flow into turbulent eddies and the formation of convective eddies are 
inseparable processes, governed by identical intertwined criteria. Marginal stability is no 
more likely to be attained in advection-dominated type flows than it is in disks. Indeed, it 
is the departure from marginal stability that is critical for maintaining vigorous turbulent 
transport, transport which is at once affected by the dynamical and the thermal properties 
of the fiow. 

The following section presents a detailed development of these results, and the final 
section of this paper discusses some the consequences for understanding the stability of hot 
plasma winds and accretion fiows. 
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2. Local Stability of a Dilute Plasma 



2.1. Physical Description 



The appearance of and T gradients as stability discriminants has fundamental 
significance. They strongly destabilize classically stable flows. Before entering a technical 
discussion, it is desirable to have a physical understanding of how this arises in some 
illustrative cases. 

In figure (1), we contrast the stability behavior of rotating fluid elements in an 
unmagnetized (la) and magnetized (lb) medium. In the unmagnetized case, the fluid 
element retains its angular momentum li on its epicyclic excursion. The angular momentum 
of an undisturbed orbit at the new location is l2- The Rayleigh criterion for stability 
is simply li < I2, so that the displaced element would drop back to a lower orbit. This 
means that the angular momentum in the disk increases outward. The magnetized medium 
is distinguished by magnetically tethered fluid elements, which enforces isorotation at 
marginal stability, i.e., Qi = Now the criterion for stability (angular momentum less 
than surroundings) is Qi < ^2- Hence, a magnetized disk requires that the angular velocity, 
not the angular momentum, increase outward for stability. This is the magnetorotational 
stability criterion. 

Figure (2) shows the conflguration appropriate to a thermally stratifled medium. In 
the unmagnetized case (2a), the fluid element retains its entropy 5*1 as it is displaced. The 
entropy of the undisturbed layer at the element's new location is 5*2. The Schwarzschild 
criterion for stability is Si < S2, which is just the condition for the blob to be denser than 
its (pressure-equilibrium) surroundings. This means the entropy increases upwards. 

In the magnetized case (2b), the displaced fluid element draws magnetic fleld lines up 
with it, and the temperature gradient now has a component along the displaced field line. 
Heat flows along this direction from the hotter to the cooler element. Under conditions 
of marginal stability, thermal conduction maintains a constant temperature between the 
fluid elements. The criterion for buoyant stability is now Ti < T2. A magnetized medium 
requires that the temperature, not the angular momentum, increase upward for stability. 
We shall now show that the classical H0iland criteria are simply but signiflcantly modifled 
to accommodate these new stability criteria. 



In the presence of a magnetic fleld, heat is restricted to flowing along lines of magnetic 
force, provided that the ion Larmor radius is much less than the coUisional mean free path 
(e.g. Braginskii 1965). This is equivalent to 



where lVci is the ion cyclotron frequency, and r is the mean free collision time. This 
inequality is generally satisfled by dilute astrophysical plasmas. Under these circumstances, 
the electron heat conduction parallel to the magnetic fleld is given by 



2.2. Analysis 




Q = -xb {b . v)n 



(2) 
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where x is the electron conductivity (Spitzer 1962), 

X ~ 6 X 10"^ Tl'"^ ergs cm"^ K"^ , (3) 

and Te is the electron temperature. & is a unit vector in the direction of the magnetic field. 
Unless otherwise explicitly stated, we shall assume that the ions and electrons have the 
same temperature T. This amounts to requiring that the mean free path A and the flow 
scale height H satisfy A/if < 1 (Cowie and McKee 1977), or 

10^^ « 1- (4) 



The usual fluid equations for a monotomic plasma are 

^ + V.(pt;) = 0, (5) 

+ {pvV)v = -V (P + - pV* + (^.v) (6) 
— ^Vx{vxB), (7) 

These are respectively mass conservation, the equation of motion, the induction equation, 
and the entropy equation. Our notation is standard: v is the fluid velocity, p the mass 
density, P the gas pressure, B the magnetic field, and $ the gravitational potential. The 
ion-dominated viscosity is ignored, as it is smaller than the conductivity by a factor of order 
the square root of the electron-to-ion mass ratio. (This assumes equal ion and electron 
temperatures.) The plasma is taken to be a perfect electrical conductor. 

We assume that in the equilibrium solution, the field fines are isotherms. Then the 
linearly perturbed heat fiux is 

6Q = -xb{6b-V)T -xb{b-V)ST. (9) 

Perturbed quantities have the WKB space-time dependence expi{k • r — cut), where k is 
the axisymmetric wavevector {kn, 0, kz)- We shall work in the Boussinesq limit, and the 
linearly perturbed equations are therefore 

V-dv = 0, (10) 

dSv , ^ 5p^^ '^^f.r. SBB\ (B-V)^„ ,,,, 
—- + SvVv^^VP--V{5P+— ^B, 11 

ot p \ 4:TT J 47rp 

= V X (Sv X B) + V X (v X SB), (12) 

ot 
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\^-^-P-8v.V\nPp-''' = ^V.8Q. (13) 
Sot p 3F 

In explicit component form, the leading order WKB equations are 

knSvR + kzSvz = 0, (14) 

-ludvR H dP-2ildv^ oT^ + 'a — 

p p^ oR A-Kp 

x{B^SB^ + BzSBz)-'^Bz6Br = 0, (15) 

- iuj6v^ + 6vr- \^ ^ + Svz R— - ik • B—^ = 0, 16 
it ait oZ Airp 

. . iA;^ 5P 6p dP ikz 

-ludvz + 2 7^^ + 1 — 

p oZ Airp 

X (5^ 6B^ + 5^ 55^) - 55^ = 0, (17) 

- icu^E^j - ifc • BSvr = 0, (18) 

dQ dQ 

- luSB^ - 6Br — — - 6Bz R— -ik-B5v^ = 0, (19) 
a In it oZ 

- itoSBz - ik • B6vz = 0, (20) 
iuj^^ + (Sv ■ V) In Pp-^/^' = ^i^ik.b {6b • V)T - (k • b)'5r) . (21) 

Note that the change in the unit vector b is 

Sb = bx(^^xby (22) 

These equations differ from the adiabatic magnetized H0iland criteria studied by 
Balbus (1995, hereafter B95) only in the final energy equation (^1]) , which contains a rather 
complicated-looking conduction term. Despite the apparent complexity, the fact that only 
the energy equation has changed allows us to simplify considerably our calculation of the 
dispersion relation, by making use of the B95 result. 

Begin by setting 6T = —T{6p/p) in the conduction term, which is an implementation 
of the Boussinesq approximation (relative changes in the pressure are much smaller than 
relative changes in the temperature or density). Then, substituting for db from equation 
(P^), and remembering that the field lines are isothermal, brings us to 

_ 2 _ 2\ Sp ^ _3^^ _ vinPp-5/3 + ^(fc • b)^.VlnT. (23) 

5 P / p 5 5P B 
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Using equations (|T8|), ( pOD and (|T^), then simplifying, leads to 

5 



5p 5vr 
P 



ILO 



where 



and G is defined in situ. 



hPuj 



5vr 39 
icu 5 



Equation (^) has two important limiting forms. When x 
the adiabatic expression 



VXnPp-^/^, (24) 



(25) 



and ^ 1, we recover 



6p 



p iuj h 



5 



(26) 



where we have introduced the radial Lagrangian fiuid displacement, 6vr = d^n/dt. This 
returns the calculation to B95. But note the result obtained if we first take the limit uj ^ 0, 
followed by X ^ 0. Then G is proportional to the ratio of the temperature-to-entropy V 
gradients, and 

^ = -^RPlnT. (27) 
P 

This implies that the u ^ limit of the dispersion relation is simply obtained by taking the 
B95 result and replacing (3/5)P In Pp~^/'^ with VlnT. Since the to ^ limit is relevant for 
flow stability, comparison of equations ( p6|) and (^) points to something remarkable: the 
stability of the flow changes discontinuously when any flnite thermal conductivity is present. 
Instead of entropy gradients serving as the discriminant for buoyant stability, temperature 
gradients are key. If the temperature — not the entropy — decreases upwards, the plasma 
will become buoyantly unstable. This is a very big difference indeed, implying that simple 
adiabatic stratiflcation (or flow) is unstable (Balbus 2000). It is wholly analogous to 
the replacement of angular momentum gradients with angular velocity gradients in a 
magnetized rotator, and the accompanying destabilization of Keplerian flow. 



2.3. Dispersion Relation 



The above reasoning suggests a simple way to obtain the desired dispersion relation: 
use 9 as a "tag," and substitute 36/5 for 3/5 where it appears as a prefactor in equation 
(2.4) in B95. This prescription gives directly 



/r2 

~4 I ~2 



36 1 

— (VP) P In + ^p(^4^2x 

5p R-^ 



where 



va = B/JA-np 



~2 2 
UJ = UJ 



-An\k • Va) 



0, 



(28) 
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We may now expand this equation using (p^), introducing a = —iuo to keep the dispersion 
coefficients real. One finds: 



T2"^ '^cond — CT 



5p V ^ 3P 



where now 



a^ = a^ + {k.VA)\ a,ond = a + ^ik-by (30) 

bP 

This is the general form of the dispersion relation. 

The marginal stability of purely evanescent modes may be studied very simply by 
passing through the point cr — > 0. The condition for stability is 

(k • va^'^ - -{VP){V\nT) - i^P(i?4^2) _ > Q ^3^) 

Note that this condition is independent of the conductivity x, for any finite value of this 
parameter. By way of contrast, when x = 0? the a —>■ result is 

(fc . vaY^ - ^(PP)(P In Pp-5/^) - ^V{R'n') - An' > 0, (32) 

which is the same as the previous condition, except for the replacement of a logarithmic 
temperature gradient by an entropy gradient. 

The condition (^) is actually much more general than the above simple derivation 
would suggest. It is in fact a necessary and sufficient condition to preclude both instability 
and overst ability. This is established in detail in the Appendix by use of the Routh-Hurwitz 
criterion. 

At this point of the calculation, the route becomes identical to B95, and the stability 
criteria are obtained by direct substitution of DlnT for (3/5)P In Pp~^/^ in equations (2.9) 
and (2.11) of that paper. For ease of reference and cross comparison, we give three forms of 
the H0iland stability criteria, and the conditions under which they are valid: 

The textbook case (Tassoul 1978) is adiabatic and unmagnetized. 

CLASSICAL H0ILAND CRITERIA. 

-l(VP).Vl„Pp-/3,i,f!||!>0. ,33) 



dP\ fOR^n? d\nPp-^l^ dR^Q'dlnPp-^/^' 
"dZJ \ dR dZ dZ dR 



> 0. (34) 



The B95 result allows for the presence of a weak magnetic field, but ignores thermal 
conduction. 
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ADIABATIC, MAGNETIZED CRITERIA. 



±(VPyV\nPp-'/-' + ^>0, (35) 
bp omR 



dP\ fdn'^dlnPp-^l'^ On'^ d In Pp-^/^' 
dz) \~dR dZ 'dZ dR 



> 0. (36) 



The result of this paper includes both the dynamics of a weak magnetic field, and the effects 
of magnetically inhibited Coulomb conductivity. 

NONADIABATIC, MAGNETIZED CRITERIA (This paper). 

/ dP\ fdQ^dlnT dQ'^dlnT\ , , 

The combined effect of a magnetic field and Coulomb conductivity is to ensure that 
free energy sources (angular velocity and temperature gradients) control flow stability. This 
result goes some way toward understanding why minimal energy and maximal entropy 
states of bound systems are associated with uniform rotation and isothermality, yet the 
classical dynamical stability criteria involve gradients of angular momentum and entropy. 
Departures from uniform rotation and isothermality are indeed a source of dynamical 
instability. It is just that magnetic tension and magnetically confined conduction are needed 
to provide the right coupling to tap into these sources. 



2.4. Some additional points 



2.4- 1- Radiative conduction 



Because of the apparent generality of the criteria (|37|) and (|38D, it is important to 
emphasize the point (Balbus 2000) that the destabilizing role of thermal conduction is 
modified substantially when diffusivity is dominated by radiative processes, as in stellar 
interiors. This form of the heat conduction is indifferent to the magnetic field. Denoting 
the radiative conductivity as Xrad, we have (Schwarzschild 1958): 



Xrad 



16TV 

Sup ' 



(39) 



where a is the Stefan-Boltzmann constant and k is the radiative opacity. When both 
Coulomb and radiative conductivity are present, the right hand side of equation ( pTD is 
modified to 

2 



3P 



ixk -biSb. V)T - {xik • bf + xradk') 6T 



(40) 
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(We assume, as before, that unperturbed field lines are isothermal. Our ultimate conclusion 
is not strongly affected by this assumption.) The dispersion relation becomes 



1^ ~4 ~2 
T2"Cr arad - CT 



—Vp(aV\nPp-^/^ + ^(k-b)^V\nT] + 
hp V 3P 



where 



2T 



rad 



0, 



(Trad = Cr+ • bf + Xradk'^)- 



The stability condition (0) becomes 



[k • va) 



e {l/p){VP){V\nT) 



l + {Xrad/x){k/k-b) 



1 - ^^(^'^^ 



4fi2 > 0. 



(41) 
(42) 

(43) 



By reducing the effective size of the potentially destabilizing temperature gradient, radiative 
conductivity is strongly stabilizing. Even when formally present, instability occurs primarily 
at long wavelengths (along the field lines) when radiative conduction is dominant, and the 
local WKB approximation we have been using breaks down. When unstable wavelengths 
are calculated to be in excess of the size of the system, the question of stability must be 
determined by a global analysis. The dynamical stability of stellar radiative interiors, 
therefore, is not threatened by this analysis (Balbus 2000). 



2.4-2. Instability in a dynamical background 

The development of a local instability in an evolving, dynamically active background 
is more complex than the analogous problem for a static equilibrium (e.g. Balbus & Soker 
1989). Both accretion flows and winds are a natural venue for the instabilities of interest, 
but we shall defer a detailed study of the dynamical linear stability theory to a separate 
investigation (Balbus 2001, in preparation). It is nevertheless both possible and useful to 
make some simple statements of a general nature. We restrict our comments to spherical 
flow. 

At t = 0, label each fluid element by a position vector r', which is a Lagrangian 
coordinate system, comoving with the unperturbed flow. The position vector r is the 
instantaneous Eulerian coordinate of a fluid element as it flows. Radial stretching of the 
flow is characterized by the scale factor 

dr 

which is a sort of Hubble parameter. The fundamental perturbation variable in a dynamical 
flow is C,r/a-, where is the radial displacement of a perturbed fluid element relative to 
the undisturbed flow. Note that we normalize C,r relative to a, since the latter tracks the 
relative separation of two points in the underlying flow. Displacements relative to a track 
true, physical changes. 
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If the growth rate of an embedded instabihty is rapid compared with the background 
flow evolution, the resulting form for ^r/o is rather intuitive: 

^ oc exp (^J\{t') dt'Y (44) 

where 7 is the value of the instantaneous growth rate obtained by ignoring background flow 
at the fluid element's position at time t'. In a spherical system, it is not difficult to show 
that the most rapidly growing thermal mode corresponds to 

^2 _ _ldPdhiT 
p dr 9r ' 

which is likely to be of order the sound crossing time. In other words, when the background 
flow is subsonic, one expects only small quantitative changes in the development of a 
perturbation from static theory. 

The breakdown starts to occur when the flow becomes sonic. Then, there is only 
one time scale in the problem, and the integral in becomes logarithmic, with power 
law behavior for the perturbations. Once the flow becomes highly supersonic, there is a 
complete breakdown of the WKB form (|4^) . Generally, such disturbances are "inflated 
away," and in any case cannot make contact with upstream fluid. Applications of the 
dispersion relation or the stability criteria to dynamically active flows should be restricted 
to their subsonic zones. 



3. Discussion 



The understanding that temperature and angular velocity gradients regulate flow 
stability in magnetized dilute plasmas has important consequences. When rotation is 
unimportant, isothermal, not adiabatic, conditions should prevail. This should be true 
whether or not field lines are free to open up. Marginal stability is a likely outcome, since 
the temperature profile is free to adjust itself while maintaining hydrostatic equilibrium. 
One interesting application is to X-ray gas in early-type galaxies. Fabian et al. (1986) 
show that the classical Schwarzschild stability criterion implies a minimum mass contained 
within a confining outer radius Tc of 

M > ^l^i^Z^, (46) 

where cq is the isothermal sound speed at observed radius tq, at the pressures Pc and Pq are 
those at Vc and tq. This constraint may be further tightened with the isothermal criterion 
to 

M > £ir^!;^. (47) 

which also has the added advantage of being rather insensitive to the pressures. (Note 
that classical cooling fiows are not affected by this instability, as the temperature increases 
outward when radiative losses are important.) 
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When application is made to a rotating systems, the notion of marginal stability is often 
inappropriate, for the simple reason that the system may not have the option of changing 
its rotation law. In Keplerian disks, the instability reduces to the magnetorotational 
instability, and that is hardly a marginally stable flow. More contemporary examples are 
ADAFs (Narayan, Mahadevan, & Quataert 1998 for a review), and CDAFs (Narayan et 
al. 2000), which were mentioned briefly in the Introduction. The latter are particularly 
striking, as proponents argue that convection itself is capable of suppressing accretion, i.e. 
that inward convective angular momentum very nearly cancels the outward transport of 
MHD turbulence. This is a rather bold extrapolation from earlier numerical simulations 
(Stone & Balbus 1996, Cabot 1996), which suggest an effective a value some three orders of 
magnitude smaller than that obtained from the magnetorotational instability. Convective 
turbulence tends to create velocity-temperature correlations that (not surprisingly) do a flne 
job of transporting heat, and a very poor job of transporting angular momentum. Recent 
MHD simulations of resistively-heated, non-radiative, quasi-spherical accretion flows (Stone 
& Pringle 2001) show no evidence of convective stalling. 

Part of the difficulty of ascribing a special role to convective angular momentum 
transport may be grasped from the form of the stability criteria appropriate to these 
flows, equations (|37|) and (|38D. CDAF proponents use, incorrectly, the classical H0iland 
criteria (|33| ) and (P4|). Among other problems, this has the effect of masking the 
dominant magnetorotational instability, which ultimately governs the nature of the angular 
momentum transport. The formal epicyclic frequency (the term proportional to the angular 
momentum gradient in ||33|) is real- valued in the flows of interest, and gives no hint of 
rotational instability. Only the adverse entropy gradients, which are clearly associated with 
convection, would seem to be involved with instability. But this is very misleading. 

The appropriate form of the instability criteria (^) and (|38| ) tell a much different 
story. The point that the thermally driven component is governed by temperature, not 
entropy, gradients is a relatively minor one in this context. (For this reason, we need not 
concern ourselves here with the possible complexities of a two-temperature plasma.) Much 
more important is the conceptual point that convective and rotational instability must be 
treated on the same footing. There is no separate "viscous" angular momentum transport 
and "convective" angular momentum transport, any more than there is turbulent transport 
and magnetic transport. There is a rotationally driven MHD instability that gives rise 
to a turbulent stress tensor, and there is nothing marginally unstable about the flow. As 
mentioned above. Stone & Pringle (2001) have performed two-dimensional calculations that 
bear this out; preliminary three-dimensional simulations are equally supportive (Hawley, 
Balbus, & Stone 2001, in preparation). 

It is, nevertheless, of great interest to understand the MHD thermal properties of 
nonrotating (or uniformly rotating) flows. What are the numerical prospects for studying 
the thermoclinic instability in nonrotating systems? It is a daunting task, because thermal 
conduction along magnetic fleld lines must be isolated, uncontaminated by numerical 
cross-fleld diffusion. When the fleld becomes highly tangled, computations become 
prohibitive. 

Even at this early stage of inquiry, something can be done. The linear instability in 
a simple vertically stratifled box with a weak horizontal fleld can be demonstrated under 
Schwarzschild stable conditions (Stone 2001, private communication). It is even possible 
to see the early stages of nonlinear instability, as shown in flgure [3]. If the effect of the 
instability is generally to comb fleld lines out, then it may well be possible to make further 
progress on the numerical front. Either by fleld line combing or by convection, the tendency 



- 13 - 



toward isothermality of a bottom-hot dilute stratified plasma seems a likely outcome. 

The results of this paper should make very clear the important conceptual point that 
even the tiniest of magnetic fields can have dramatic consequences for the macroscopic 
stability of astrophysical plasmas. Classical hydrodynamical results can be qualitatively 
incorrect, and great care must be taken before uncritically taking them over into magnetized 
systems. Subtlety need not imply great complexity: the relative simplicity of our 



fundamental results, equations (|37|) and (^), is encouragement that the most important 
MHD stability properties can also be conceptually simple. And since it seems certain that 
the final word has not been said on the topic, more surprises are likely. 

It is a pleasure to thank J. Hawley and C. Terquem for detailed comments on a 
preliminary draft of this paper which led to significant improvements in the presentation, 
and M. Fall for pointing out to me that the new stability criteria sharpen dark matter 
constraints in X-ray clusters. I am also most grateful to J. Stone for providing me with 
fig. [3], and for discussions related to the numerics of magnetized thermal conduction. This 
research was supported by NASA grants NAG 5-7500, NAG5-106555, NAG5-9266, and 
NSF grant ASTOO-70979. 



Appendix: Necessary and Sufficient Criterion for Stability 

To prove that the condition ( pTD is a necessary and sufficient condition for any type of 
instability, we first expand the polynomial in equation (29) in the form 



aoa^ + aicr^ -|- a2(T^ + a^a"^ + a^a + = 0. 



(48) 



where 



ai 



2xT P 



5P kl 



k • by 



as 



02 

2xT 



k'z 



3 



2(fc • va?^ - ^VPVlnPp-^/^ - ^V{R^n'' 

k^ 1 „ 1 



05 



5P 

a4 = (fc • vaY 
2xT 



k . by 



2(fc • va)''^ - -VPVliiT - -^V{R*n^] 
kz P ^ 

\2'^ '-' T-. o T-. i„ o „-5/3 'T-.^D4r,2^ 



{k ■ vaYtj - —VPVlnPp-^'' - —ViR'n') 

K r7 



hp 



5P 



k • by{k • VAy 



{k • VAj 2 



k'^ 1 

-tVPVhiT 



P 



1 

R^ 



in' 



To make things more compact and manageable, notice that 



(49) 
(50) 

(51) 

(52) 

(53) 

(54) 



.ApppinPp-5/3 + --DPVlnT = -^{VPf > 0, 
bp p 5pP 



(55) 
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and if we furthermore define 
then we may write 



(56) 



as 



5P 



[k ■ by{a2 



{k • VAfia2 - 6), 05 



2xT 
5P 



{k-bfik-VAfia2-S-e). 

(57) 

For stabihty, we require that the real part of a must be less than zero, i.e., that the 
roots of equation ( ^Sf ) must all lie in the left complex a plane. Polynomials with this 
property are known as Hurwitz polynomials (e.g. Levinson & Redheffer 1970). If the left 
side of equation ( ^8]) is a Hurwitz polynomial, both pure exponential instability as well as 
overstability are precluded. This will be so if and only if 





ao 











as 


a2 


ai 


ao 





as 


a4 


as 


a2 


ai 








as 


a4 


as 














as 



> 0, 



(58) 



and all five determinants obtained by expanding along the diagonal, i.e.. 



ai. 



(l\ (Xq 

as a2 



etc. 



(59) 



are each individually > (Levinson & Redheffer 1970). We shall denote the determinant of 
each n X n matrix thus obtained as det(?T,), so that 

det(l) = ai, det(2) = aia2 — Ooas, 

and so forth. The requirement that the left side of equation (^) be a Hurwitz polynomial 
may then be succinctly stated as 



det(n) > 0, 



n 



1,2. ..5. 



(60) 



This is the Routh-Hurwitz criterion. The criterion (U) in the text amounts to the condition 
as > 0, which follows directly from the two demands det (4) > and det(5) = as x det(4) > 0. 
We shall now prove that the general Routh-Hurwitz criterion is satisfied for our dispersion 
relation. 



The first non-trivial requirement is 



det(2) = aia2 — aoas > 0. 



(61) 
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The determinant is easily calculated, 

det(2) 



\kzj 5P 
and is, indeed, positive. 

The calculation for det(3) is a bit more complicated. 



(62) 



det(3) = as det(2) — ai 



CI5 CI4 



(63) 



which simplifies to 



det(3) 



^kzj \5PJ 
But the term in square brackets is just 



— -r-fr- {k • b) e a2 - e - [k • va 



(a2 - 5 - e) + 4fi' + [k ■ va^ tt " 1 • 



(64) 



(65) 



According to equation (0), the sum of the first three grouped terms must be positive if 
as > 0, and the remaining two terms are each manifestly positive. Hence det(3) > 0. 

The calculation for det(4) proceeds as follows: 





ai 


ao 





det(4) = a4 det(3) — 


as 


a2 


ao 




as 


a4 


a2 



(66) 



The determinant cofactor of as in the above may be expanded and simplified, reducing to 





ai 


ao 







as 


a2 


ao 




as 


a4 


a2 


Using this in equation 


m 


for 



(67) 



det(4) 



kz, 



2,2 



{k ■ b)'' {k ■ VAYe 



(68) 



which is manifestly positive. Finally, as noted, det(5) = as det(4) > 0, if as > 0. We 
have therefore shown that the condition as > 0, equation (Rip in the text, ensures that 
the real parts of the roots of the polynomial in equation (29yall lie in the left half of the 
complex a plane, and that the criteria (B7[ ) and (|55| ) are necessary and sufficient for the 
convective- rotational stability of a hot plasma. QED. 
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(lij h 

I J < 1 2 for Stability 




UNMAGNETIZED 




MAGNETIZED 



Fig. 1. — (a). In an unmagnetized disk, specific angular momentum is conserved. Stability 
may be determined by comparing the angular momentum of an outwardly displaced fluid 
element (Zi) with that of its surroundings {I2). If h < h, the displaced fluid element has 
less angular momentum than it needs to remain in its new orbit, and it drops back. This is 
equivalent to the Rayleigh condition dl/dR > 0. (b). In a magnetized disk, it is the angular 
velocity of tethered fluid elements that tends to be conserved near marginal stability. Hence, 
if Qi < Q2, the displaced fluid element has less angular momentum than it needs to remain 
in its new orbit. This corresponds to the MRI stability condition, dfl/dR > 0. 
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Sj < S2 for stability 




UNMAGNETIZED 



b) 






1 rrt n-2/5\ rwi -^-2/5 




/ TjKT^ for stability 


It p'^^^ I 

v J 






MAGNETIZED 



Fig. 2. — A process similar to that depicted in fig. [1] occurs in a thermally stratified 
disk. (a). In an unmagnetized layer, entropy is conserved. Stability may be determined 
by comparing the entropy of an outwardly displaced fluid element {Si) with that of its 
surroundings (5'2). If S-i < 82, the displaced fluid element has less entropy than it needs to 
remain in its new orbit, and it drops back. This is equivalent to the Schwarzschild condition 
of increasing upward entropy, (b). In a magnetized conducting layer, it is the temperature 
of tethered fluid elements that tends to be unchanged near marginal stability. As in (a), 
if Si < S2, the displaced fluid element has less entropy than it needs to remain in its new 
orbit. But since S oc In TP-2/5, and pressure balance is maintained, this corresponds to a 
thermoclinic stability condition: the background temperature must increase upwards. 



-19- 




Fig. 3. — Development of thermoclinic instability in a Schwarzschild-stable layer. Magnetic 
lines of force are shown after one Alfven crossing time, initial seeding with rms 1% random 
initial vertical velocity perturbations. Initial thermal energy density is 1600 times magnetic; 
initial field lines are isothermal and horizontal; vertical grid runs from z = 1 to 2, 
initial temperature profile is gravitational field is l/z^] x is 0.05; grid is 128 x 64. 
The interpenetration of cool and warm blobs appears very similar to classical convective 
instability. 



